Transport properties of individual C 60 -molecules 
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Electrical and thermal transport properties of Ceo molecules are investigated with density- 
functional-theory based calculations. These calculations suggest that the optimum contact geometry 
for an electrode terminated with a single- Au atom is through binding to one or two C-atoms of Ceo 
with a tendency to promote the sp 2 -hybridization into an sp 3 -type one. Transport in these junctions 
is primarily through an unoccupied molecular orbital that is partly hybridized with the Au, which 
results in splitting the degeneracy of the lowest unoccupied molecular orbital triplet. The trans- 
mission through these junctions, however, cannot be modeled by a single Lorentzian resonance, as 
our results show evidence of quantum interference between an occupied and an unoccupied orbital. 
The interference results in a suppression of conductance around the Fermi energy. Our numerical 
findings are readily analyzed analytically within a simple two-level model. 



PACS numbers: 73.63.Rt, 73.63.-b, 73.23.Ad 

I. INTRODUCTION. 

Ceo on metal surfaces is an important model system for 
understanding basic processes in binding of (conjugated) 
organic molecules to metal electrodes, and has been stud- 
ied in the past both experimentally and theoretically?" 8 -! 
Transport through Ceo has motivated investigations by 
experimentalists and theorists from early on with an em- 
phasis on correlation physics like the Kondo effectP^ or 
vibrational degrees of freedom^"^. 

Concerning the linear conductance, the situation 
seems to be particularly well investigated with Cu- 
electrodes, where a combination of ab-initio based cal- 
culations and STM-experiments have provided a detailed 
under standing P^EH These investigations suggest that the 
conductance, G, of Ceo bound to Cu(lll)-substrates is 
sensitive to the anchoring mechanism. In general, G 
is relatively large for a single Cu-atom contacting Ceo 
immobilized on a Cu(lll) surface (~0.13 Go, where 
Go = 2e 2 /h). 22 In contrast, the conductance of Ceo 
on Au-electrodes has not been studied extensively. Me- 
chanical break junction experiments have found that the 
conductance of a junction can be as high as 0.1 — 0.2 
q J23 | 24 | w \ l [\ e STM-break junction experiments report 
much smaller values with a very broad scattering in con- 
ductance histograms for Au-, Pt- and Ag-electrodes 25 . 

Understanding the transport mechanism in Ceo/Au 
junctions is of interest for two reasons. First, any anal- 
ysis and design of transport processes through single 
molecules relies upon an understanding of the influence 
of the electrodes! 2 ^ Second, the specific molecule Ceo- 
molecule might play a special role in the context of con- 
tact formation because it was recently proposed to be a 
suitable general a nchor group due to its size and elec- 
tronic conjugation! 27 * 28 * 



In this paper we address the transport characteristics 
of single- molecule junctions formed using Ceo molecules 
attached to Au electrodes. Our calculations are em- 
ploying the DFT-based non-equilibrium Green's func- 
tions (NEGF) formalism^. Our first theoretical result is 
that Au-electrodes are invasive: Au single adatoms have 
a tendency to form chemical bonds with C-atoms and 
thus locally affect the sp 2 -conjugation in Ceo, similar to 
what has been reported for Cu electrodes. We find that 
the alternative scenario, where the adatom resides in a 
hexagonal/pentagonal facet of Ceo is not energetically fa- 
vored as it has a binding energy that is 0.5eV lower. We 
thus confirm statements reported in Ref. [30j Our second 
result is that Ceo has a slight tendency to charge neg- 
atively on Au with single adatom binding, though not 
as much as on Cu (or Ag)P^ This contrasts results on 
flat Au(lll) where the charge transfer is negligible!^ 
As a consequence, with Au-electrodes transport is more 
LUMO-dominated. However, in contrast with Cu, it 
will in general not be close to resonant. Third, we 
determine the transport characteristics of several Ceo- 
j unctions that exhibit somewhat different contact ge- 
ometries by calculating the transmission function, T(E), 
and the Seebeck-coemcients, S(E). For two geometries 
data similar to ours have been reported before 32 . Our 
analysis goes significantly beyond earlier work because 
we provide evidence that interference between transport 
channels plays a quantitatively important role, especially 
when the Fermi-energy is situated between HOMO- and 
LUMO-resonances of the Ceo molecule. We show that 
deviating from earlier claims 33 the conductance at the 
Fermi energy is not resonant but rather strongly sup- 
pressed due to destructive interference from two strongly 
coupled transport channels. This suppression leads to a 
sharp, step-like increase of the Seebeck-coemcient near 
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FIG. 1: Energy profiles from DFT with Grimme corrections: Adatom geometries with Cu-electrodes (a) and Au-electrodes (b). 
(c): Flat Au-electrodes facing Ceo hexagons. 



the minimum transmission energy. We analyze our find- 
ings within an analytic model for a two-level system and 
find that it supports the results of DFT-based transport 
calculations. In addition, this analysis also shows that 
due to cancellation effects, electrical currents driven by 
heat gradients (rather than by a bias voltage) remain 
almost unaffected by destructive interference effects. 



constraints, adatoms may also sit on-top of C-atoms or 
in penta-hexa bridge position. Since the energy differ- 
ence between these three geometries is under 100 meV 
even close to the minima, these three geometries should 
be energetically accessible, especially in break-junction 
experiments where the electrode structures can deviate 
from pyramidal tips. 



II. SYMMETRIC CONTACTS WITH SINGLE 
ADATOM 

A. Method 

We compute the total energy of the extended molecule 
(Ceo plus contact model) for different contact geome- 
tries using density functional theory (DFT). To this end, 
we employ the TURBOMOLE package^ with the BP86 
functional 35 and van der Waals interactions included on 
the level of Grimme corrections^ see appendix [a] for 
details on them. The geometry was optimized in the fol- 
lowing way: the relative position of all electrode atoms 
was fixed with bond-lengths as given by Au-bulk. The ge- 
ometry of Ceo was optimized in vacuum. Thereafter, Ceo 
was approached to the electrode with different surface- 
molecule contact geometeries: hexagon or pentagon on 
the flat surface or atop (C-Au-atom) and bridge position 
for single-adatom geomeries, see Fig. |TJ For each contact 
geometry a trace binding-energy was recorded as a func- 
tion of distance. In this process the geometry of Ceo was 
fixed. 



B. Results: Adatom geometries 

Binding energies for different inter-electrode distances 
are displayed in Fig. [TJ One infers from Fig. |l^L,b 
that Cu or Au-adatoms prefer bonding in bridge posi- 
tions, where two hexagons touch each other (hexa-hexa- 
bridges). However, for larger distances between the elec- 
trodes, i.e. more generally for relaxation under external 



Discussion: Carbon atoms of Ceo in vacuum are sp 2 - 
hybridized. However, the bond angle (108°) is relatively 
far from the "ideal" sp 2 -hybridization value (~ 120°). 
Therefore, it is plausible that Ceo is susceptible for bond- 
ing with adatoms 37 , in contrast to simple sp 2 -hybridized 
carbon like graphene or other conjugate molecules. Fig. 
|l^,b reveal that the binding energy of Cu to Ceo exceeds 
that of Au by 0.5eV. On a qualitative level this observa- 
tion could be related to the fact that the KS-work func- 
tion of Cu is situated considerably above the one of Ceo 
in vacuum (see Fig. |2|. Therefore, a moderate flow of 
electrons into the Ceo-LTJMO might increase the inter- 
action with the Cu-surface as compared to the one of 
Au. 



The relatively strong interaction between the carbon tt- 
system and the Au-adatom also manifests itself in the ra- 
tio of the GGA binding energy (ignoring Grimme the cor- 
rections) to the total binding energy, p = Eqqa/ (Eqq\ + 
^Grimme)- It helps to quantify how close the bond is to 
being covalent. If p « 1, the van der Waals contribution 
is negligible and a covalent chemical bond has formed. In 
contrast, if p <C 1, bonding is predominantly of the van 
der Waals type and the adatom should be thought about 
as being physisorbed. 

The ratios p are given in Table [I] along with the cor- 
responding total binding energies and bond lengths. For 
both Cu and Au-electrodes, the values p > 1/2 suggest 
that bonding is predominantly covalent consistent with 
the expectations formulated in the preceeding paragraph. 
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Electrode 


Position 


^bond [A] 


£bind [eV] 


P 


Cu 


C-atom 


2.0 


-1.83 


66% 


p-h bridge 


2.2 


-1.89 


61% 


h-h bridge 


2.2 


-2.05 


63% 


Au 


C-atom 


2.2 


-1.29 


56% 


p-h bridge 


2.4 


-1.31 


54% 


h-h bridge 


2.2 


-1.49 


54% 


Au (flat) 


hexagon 


3.0 


-1.02 


-32% 



TABLE I: Characteristics of lowest energy molecular junction 
geometries where /bond is the Au-C or Cu-C bond length, 
^bind the binding energy per electrode and p = i^GGA/^bind 
gives an indication of the covalent contribution to the bond. 



-2.5 
-3 

-3.5 
-4 

£ -4.5 



-5 



-6.5 
-7 
-7.5 



LUMO+1 (3X) 



LUMO (3X) 



HOMO (5X) 



HOMO-1 (4X) 



C. Results: Flat Au(lll) electrodes 

We also consider flat Au(lll) electrodes (without an 
adatom) facing a hexagon of Ceo- The binding energy 
profile of that configuration is displayed in Fig. [J. The 
bond distance here is measured from the position of the 
nuclei of the first Au(lll) layer. Most notably, the pure 
GGA functional without Grimme corrections would in- 
dicate that the configuration at a distance of ~ 3.0 A is 
non-binding. This is reflected in the sign of p in Table [T[ 
The bond has no covalent contribution in this geometry, 
consistent with the higher coordination number of sur- 
face atoms compared to adatoms. As shown in Table 
|lj the bond distances for the C6o-Au(lll) geometry is 
0.8 — 1.0 A larger than in the case of Ceo-adatom geom- 
etry. The binding energy is of similar magnitude as in 
the adatom case, indicating that this configuration may 
indeed be relevant in experiments. 



FIG. 2: Kohn-Sham energy levels of Ceo in vacuum (black 
lines). Fermi energies of Cu (red) and Au (yellow) estimated 
from DFT calculation for the 24-atoms clusters used for the 
transport calculations. 



a discussion of this discrepancy, see Sec. IV E| ) The po- 
sition of the chemical potentials of Au and Cu relative 
to Ceo energy levels are also shown in Fig. [2] These 
values are from DFT calculations on 24-atoms clusters 
since the same values will be used in our transport calcu- 
lations. The position of these levels suggest that charge 
transfer between the Au-electrode and molecule will be 
relatively week, in comparison with Cu-electrodes where 
the molecule can pick up a pronounced negative excess 
charge. 



III. ELECTRONIC STRUCTURE 

The nature of the molecule-electrode bond is of cru- 
cial importance for Ceo transport properties. Since the 
formation of a covalent bond with Ceo implies that the 
conjugation of the 7r-electron system is broken at the con- 
tact, a transport barrier forms and the molecule should 
be considered as "weakly coupled" . In this section we dis- 
cuss the effect of the bonding on the electronic structure 
of the molecule. 



A. Free molecule 

The position of relevant molecular orbitals of Ceo in 
vacuum, as obtained from DFT calculations, is shown 
in Fig. [2] The Highest Occupied Molecular Orbital 
(HOMO) and Lowest Unoccupied Molecular Orbital 
(LUMO) levels are 5-fold and 3-fold degenerate. The 
associated HOMO-LUMO gap is A - 1.6 eV, slightly 
underestimating the experimental value 2.3 eV^. (For 



B. Local Density of States 

When the molecule is in contact with metal electrodes, 
the position of these molecular levels will shift and ex- 
perience a lifetime broadening. The electronic structure 
of the junction is represented by the density of states 
projected on the Ceo (local density of states, LDOS). 
We calculate it for the geometries of lowest total energy 
with the DFT-based Green's function formalism for non- 
interacting particles described in Appendix [B] and Ref. 
US The result is displayed in Fig. [3^,b. 

To quantify further, we parametrize the LDOS as a 
sum of Lorentzians, 

Values of the fitting parameters, resonance position e n 
and broadening 5 n , are given in Table |H| for junctions 
where the adatom is facing a C-atom of the Ceo (Fig. [i]) . 
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FIG. 3: (Color) Local Density of States (LDoS) projected on Ceo for Cu-electrodes (a) and Au-electrodes (b) for the three 
geometries of lowest total energies identified in Fig. [I] (c) LDoS of Ceo in a junction with direct coupling to a flat Au(lll) 
surface for the optimum geometry, see Fig. [I]:. 







HOMO 


LUMO 


LUMO+1 


Au 


e n - E F (eV) 


-1.68 


-1.38 


-1.37 


-1.35 


-1.35 


0.10 


0.25 


0.28 


1.15 


1.28 


1.33 


S n (meV) 


149 


26 


8 


1 


6 


116 


8 


6 


175 


21 


1 


Cu 


e n - E F (eV) 


-1.59 


-1.48 


-1.48 


-1.45 


-1.44 


0.06 


0.12 


0.16 


1.02 


1.15 


1.24 


S n (meV) 


55 


15 


17 


5 


3.5 


50 


12 


11 


60 


40 


1.1 


Au (flat) 


e n - E F (eV) 


-0.78 


-0.77 


-0.76 


-0.75 


-0.75 


0.84 


0.86 


0.87 


1.92 


1.92 


1.93 


S n (meV) 


114 


116 


4 


2 


2 


82 


75 


17 


9 


42 


53 



TABLE II: Linewidths, S n , and positions of energy levels, e n , for the geometry where the adatom sits on-top of a C-atom, and 
the flat surface geometry. Levels are classified according to their position in the energy spectrum of the molecule in vacuum 
(HOMO, LUMO, HOMO+1). 



1. Adatom geometry 

The splitting of molecular energy levels in Fig. |3j es- 
pecially the LUMO ones, indicates that the formation 
of the chemical bond constitutes a significant perturba- 
tion in the sense that one of the LUMO states of the 
molecule splits away from the others. It hybridizes more 
strongly with the electrode states as exhibited by the cor- 
responding increased level broadening. Electron trans- 
port through Ceo will be mostly via this level. Table [TT] 
reveals the effect of electrode coupling on level splitting 
and broadening is about two times stronger for Au- than 
for Cu-electrodes. 



2. Flat Au(lll) contacts 

With a flat Au(lll) surface (see Fig. [J) £Vermi resides 
in the middle of the HOMO-LUMO gap, much closer to 
its vacuum position (see Fig. [2| than in the adatom ge- 
ometries. This indicates, that partial charge transfer is 
weaker with a flat electrode as compared to the case of 
adatoms, reflecting that binding is purely of the van-der- 
Waals type. Consequently, the splitting of the energy lev- 
els (see Table |TT|) is much smaller than that of the adatom 



geometry. In other words, the presence of the electrodes 
implies only a weak symmetry breaking with an overall 
small effect on the molecular frontiers orbitals. We as- 
sociate the large quantitative differences in the observed 
shifts and broadenings with wavefunction overlaps: sym- 
metry related extinction results in reduced hybridization 
matrix elements. 



IV. TRANSPORT CALCULATIONS FOR C 60 

If we were to consider a system with N channels that 
do not interfere with each other, then by definition the 
transmission (per spin) could be written as a sum of 
Lorentzians, 

with parameters given in Table [TTJ The formula yields a 
good approximation, usually, in the presence of symmet- 
ric coupling if a single transport resonance dominates. 
However, in cases where transfer amplitudes of different 
channels can be of comparable magnitude interference 
terms may become significant and the approximation pi) 
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breaks down. Then, in principle, the Landauer formula, 



Eq. (B4) in the appendix, has to be employed. 



A. Results: Transmission function 

The full transmission function, T(E), as obtained from 
the Landauer formula in the NEGF-formulation 29 is dis- 
played in Fig. [4^i, b for the Au-C atop-geometry. It ex- 
hibits pronounced non-Lorentzian features, with a sup- 
pression of the transmission in the valley region. Com- 
paring with the model of isolated resonances, Eq. ([2]), one 
can see that the missing cross-terms between different 
transfer modes explains this behavior. We conclude that 
interference between transport resonances play a quan- 
titatively important role in electron transport through 
Ceo when the chemisorbed Ceo-molecule is only weakly 
charged. The conductance could be reduced by roughly 
one order of magnitude as a result 



B. Discussion: effective two level model 

We interpret our findings for the transmission function 
using an effective two level model. Its precise definition 
together with a deriv ation of basic properties are given 
in the appendix |Dp^ The salient feature of the simplified 
model are summarized by the following set of equations: 



T(E) = T 



Ti±T i 

(E-e )(E-e 1 ) 



•7o7i 



V(£-eo) 2 +7oV(^i) 2 +7i 2 



(3) 
(4) 



The two levels that one should refer to here derive from 
the C 60 -HOMO quintett and LUMO-triplett, that ex- 
hibit resonances with the strongest broadening. Accord- 
ingly, we see from Table |TT] the set of parameters: eo = 
-1.68eV, ei = O.leV, 70 = 0.149meV, 71 = 0.116meV. A 
decomposition of the LDoS given in Appendix [C] further 
substatiates this simplification. 

The sign in Eq. ([3| controls the effective mixing be- 
tween the two transport channels, i.e. whether they in- 
terfere constructively (minus sign) or destructively (plus 
sign) energies in the valley region, eo < E < e\. As we 
explain in the appendix, destructive interference occurs 
in two-level models where both states couple with simi- 
lar strength to both reservoirs. For Ceo, we expect that 
there should not be any important difference between the 
coupling of the HOMO and LUMO levels to the leads, we 
can expect to see destructive interference. 

Indeed, already from Fig. [4]we can see that for the case 
of the C6o~Au-j unction, the transmission in the valley 
region is very strongly suppressed supporting our claim. 

Hence, we conclude that the plus-sign should be chosen 
in Eq. (|3|. Furthermore, with symmetric coupling we also 
have 



Ti(E) 



It 



v2 ' 



0,1. 



(5) 



The sign in Eq. ([3| is the toy model's only ingredient in 
T(E) that is not fixed by the LDOS alone. 

In Fig. [5] (left) we compare the full transmission with 
the one from the toy model, Eq. Indeed, in the 

valley region the toy model reproduces the transmission 
and all its non-Lorentzian features well up to a small 
shift of the minimum-transmission energy. This shift 
is readily explained, e.g., by residual energy dependen- 
cies in the pole-positions due to the structured density 
of states in the reservoirs. For the minimum conduc- 
tance at energy e* we obtain a parametrical estimate 
^constr ~ 47071/ 1 eo — e*||ei — e* | for constructive interfer- 
ence and 



destr 



1 9 

4 co 



(6) 



in the other case. 



C. Thermopower 

We complete our account of Ceo transport proper- 
ties with a discussion of the thermopower of an Au-C6o- 
Au junction. The thermopower or Seebeck coefficient S 
determines the magnitude of the built-in potential de- 
veloped across the junction when a temperature differ- 
ence AT is applied.^ In particular, molecular junction 
thermopower can be useful in determining the dominant 
molecular orbital for transp ort and the identity of the pri- 
mary charge carriers ! 4Q I 41 1 With the additional pre sence 
of an external voltage bias AV across the junction^ 42 * 43 *, 
the total current / in this case is simply 



2e 2 

T [-T(^ F )A^- 



■T(E F )S(E F )A%] (7) 



where T(E-p) denotes the transmission at the Fermi 
energy. For non-interacting electrons, the Seebeck- 
coefficient, S(E), is closely related to the transmission 
function T(E), ultimately, because all transport pro- 
cesses are controlled by the tunneling probabilities of 
electrons with a given energy through the barrier. We 
have (using the convention e = |e|) 



S(E) 



Tr 2 kl1d\nT{E) 



3e 



dE 



(8) 



where T denotes the (electronic) temperature. The sys- 
tem specific information is all encoded in the logarthmic 
derivative, which we now discuss. 

Fig. [SJd displays the logarithmic derivative for our 
model system, Ceo- Again, we can convince ourselves 
that the two-level model with destructive interference ac- 
counts well for the salient features of the full DFT-based 
trace. A striking characteristics to be observed here is the 
step-like transition, changing sign, that the derivative un- 
dergoes when the energy sweeps by the valley minimum 
point e* ~ — 0.51eV. We can estimate the parametrical 
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E-E F [8V] E-E F [eV] E-E F [eV] 



(a) (b) (c) 

FIG. 4: (Color) Transmission functions for different electrode models. Adatom in atop-position with Cu-electrode (a), and 
Au-electrode (b). (c) diplays the flat Au(lll)-surface. Red line: DFT-based transport code with Landauer formula/NEGF- 
formalism. Green line: isolated resonances model, Eq. Q. Plot highlights the effect of crosstalk between different transport 
channels in the transmission valley regime of energies. 




E-E F [eV] E-E F [eV] E-E p [eV] 

(a) (b) (c) 

FIG. 5: (Color) (a) Comparison of transmissions obtained from the transport code (Fig. J2Jb), solid red) and the two-states 
model (solid blue), see Eq. |3|. In order to highlight the impact of the interference term, T i, Eq. |3|, results for constructive 
and destructive situations are given, (b) Logarithm of the transmission function (dashed lines) and its derivative. The latter 
represents the system specific information content of the Seebeck-coefficient, Eq. (c) System specific, energy dependent 
characteristics, T(E)d\nT(E)/dE, of the thermal current Eq. Plot highlights a result of the two- level model: traces for 
destructive (blue) and constructive (brown) interference give nearly coinciding results, even though their electronic transmission 
deviates by orders of magnitude, see the transmission functions given in the left figure. 



dependency of the slope at e* employing Eqs. (|3]|5|): 



dlnT(E) 



dE 



(£- € *)( €l - €o ) 
27071 |A Ai| 



\E ~ e*| < 7o, 7i 
(9) 

where A$ = — e* . For the typical cases where 70 and 71 
are comparable, e* is halfway between e\ and eo, so the 
expression simplifies further: d\nT(E)/dE\ e * « 2(E — 
e*)/7o7i. Hence, the slope diverges in the weak coupling 
limit, where 70 1 tend to zero at fixed level splitting |ei — 
eol- ' 

The estimate for the slope is valid for energies in a 
vicinity of width 70,1 about e*. Hence, the logarithmic 
derivative takes very large magnitudes ~ ±1/70,1 near 
the center of the valley region. It is only near the reso- 
nances where it reaches similar values, e.g. ~ I/71 near 
e\. The intermediate region interpolates between these 
to maxima. 



This behavior is typical of systems that exhibit almost 
perfect destructive interference, so that T(E) approaches 
zero near some energy e*. It is completely absent with 
constructive interference, see Fig. [5] (center). There the 
logarithmic derivative has a parametrically small slope of 
the order of ~ l/|AoAi| (rather than I/7071) near e*. 

The temperature driven current is, up to system un- 
specific prefactors, given by the product Td\n(T)/dE. 
Each factor has been seen to be very sensitive to the 
sign m Eq. {§) in our discussion. For the product this is 
not the case, as one infers from Fig. The reason 

is that the suppression of T* for the case of destruc- 
tive interference is largely compensated by the strong 
slope in the logartithmic derivative. Indeed, in the 
valley region we have for the product with either de- 
structive of constructive interference a similar behavior: 
T*dln(T)/dE « (e - e*) 7 2 /J 4 , where 7 = 7o = 71, £ = 
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A = — Ai = (eg — ei)/2 has been assumed. 

Hence, we arrive at the following conclusion: In 
the presence of destructive interference (minus sign in 
Eq. (|3|) the bias voltage driven current can be suppressed 
by orders of magnitude in the valley region. Neverthe- 
less, an electrical current driven by a thermal bias reaches 
similar values as it would in the absence of interference 
effects. 



D. Contact geometries revisited: asymmetry and 
chain formation 

In this section we show, that the shape of the transmis- 
sion function can be significantly modified by changing 
the contact geometry. 



T 




E-E F [eV] 



1. Single adatom contacts breaking inversion symmetry 

By selecting different pairs of C-atoms on the molecule, 
one modifies the phase difference between parallel trans- 
mission paths (Fig. [6|. The situation is similar to the 
case considered perviousl y 44 1 45 1 , except that there a torus 
geometry was considered while we investigate a sphere. 
The different contact geometry has an impact on the ef- 
fect of Toi. Specifically, the transmission function for 
the second geometry (maroon trace in Fig. [6| is very 
similar to a pure superposition of Lorentzian resonances, 
suggesting that the Tqi in Eq. [3] does not contribute sig- 
nificantly and interference effects between parallel paths 
are much weaker than in the other cases. 

In contrast, destructive interference reappears in ge- 
ometry (3), Fig. [6j This geometry is however, far from 
perfectly symmetric, and therefore the simplified model 
Eq. Q cannot be expected to hold. Instead, in prin- 
ciple the m ore complete formula given in the appendix, 
( D13|D15| ), should be applied. Indeed, the fit based on 
the simplified expression (3]) does not properly reproduce 
the behavior of T(E) in the valley region. (Based on the 
discrepancy one expects cos(\I/ c ) ~ —0.8.) 

We emphasize that there is a very large variability 
of the conductance in the valley region even though all 
electrode positions (1-3) are associated with similar res- 
onance positions and broadenings. We attribute this be- 
haviour to the mixing angle \I/ C oscillating from ir (ge- 
ometry 1) to (2) back to larger values \I/ C « 2.5 (3) 
again. This observation we take as support for our claim 
that the variations observed in the transmission function 
are due to a modification of the phase difference between 
parallel paths. 



2. Au-contact chains 

In order to investigate the development of interference 
with decreasing level broadening, we consider here ge- 
ometries where the molecule is included between Au- 



FIG. 6: Transmission function for different anchoring ge- 
ometries: (1) Symmetrical geometry from previous plot Fig. 
[5] Asymmetrical geometries (2, maroon; 3, blue). Posi- 
tion (2) can be fitted by adding two Lorenzians (dashed line, 
70 = 0.106,eo = —1.6,71 = 0.096,ei = 0.13), indicating that 
interference effects are weak. In contrast, geometry (3) ex- 
hibits destructive interference. The model (|3| fails in the 
valley region (blue dot dashed line) , mainly because the angle 
^ c (defined in Appendix [dJ is not close to tt. From the blue 
data trace one estimates roughly cos ^ c « —0.8 together with 
an asymmetry 7ol/7or ~ 0.45. 



chains. In this configuration, the number of incoming 
and outgoing lead channels is limited to essentially a sin- 
gle one. This reflects in the local density of states at 
the chain terminating Au-atom that the molecule cou- 
ples to. It is more strongly structured as compared to 
the case with a single Au-adatom, only; the number of 
states that are ready to hybridize with the C6o~orbitals is 
reduced. As a consequence, broadening of molecular or- 
bitals contacting Au-chains is in general weaker, and also 
more complicated since a convolution of two structured 
functions (LDoS on molecule and contact-atom/ Au- wire) 
is involved. 

The transmission functions obtained for one, two and 
three-atoms chains geometries shown in Fig. [7] support 
these expectations. We observe a progressive develop- 
ment of large amplitude anti-resonances with increasing 
chains length. They reflect the fact molecular states and 
wire states can cooperate in a complicated manner which 
allows, in particular, for more levels to develop interfer- 
ence patterns in valley regions. In this way, transmission 
values below 10 -4 can come about for 3-atoms chains, 
which suggests that extremely low conductance values, 
between 10 _3 Go and 10 _4 Go, see Fig. [7| can be ob- 
served with this type of junctions. This finding becomes 
particularly interesting in view of the fact, that a sin- 
gle Au-chain is well known to exhibit a single perfectly 
transmitting channel. Our result Fig. [7] gives a solid 
demonstration, that due to quantum effects even a per- 
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symmetry and the fact, that the HOMO- and LUMO- 
orbitals are conjugated, coupling well in similar ways to 
both leads. 

Due to the large, interference induced slope of T(E) in 
the valley region, one could expect to observe very large 
conductance fluctuations in the experiments - despite of 
the high molecular symmetry of Ceo ~ due to weak en- 
vironmental capacitive couplings. Our study would also 
suggest that the Seebeck coefficient does not exhibit such 
strong fluctuations because its dependency on the level 
alignment in the valley region is relatively weak, see e.g. 
Fig. [5] center. Both qualitative features are indeed ob- 
served in Ref. |25j On a quantitative level, we observe 
that the theoretical estimate for the Seebeck-coefflcient 
overshoots the experimental one roughly by a factor of 
2-4. 



FIG. 7: Transmission functions for contacts made via adatoms 
(1), two- atoms Au-chains (2) and three- atoms Au-chains (3). 



V. CONCLUSIONS 



feet conductor can be a very invasive means to fascilitate 
an electrode coupling. 



E. Additional remarks 

We add several remarks on artifacts of DFT-based 
transport calculations and about experiments. 

First, the functionals used in our study are well 
known to underestimate the true HOMO-LUMO gap, 



as was pointed out already in section |III A| The miss- 
ing derivative discontinuity influences the alignment of 
the molecule-based and metal-based electronic levels; in 
general it leads to an overestimation of charge transfer. 
Hence, while one expects that the qualitative features of 
the DFT-based transport calculations are captured cor- 
rectly, the positioning of the LUMO with respect to the 
Fermi-energy, Ep, should be slightly too close. The true 
Ey is probably situated somewhat closer to the valley 
region than seen in the DFT-calculation. 

Second, recent research has shown that DFT-based 
transport calculations employing exact functionals repro- 
duce the exact transmission for interacting single level 
models.^E^ At present, a rigorous generalization of the 
statement to models with several levels does not exist. 
In constrast, the numerical results of Ref.™ show devia- 
tions between exact conductances and DFT-based trans- 
port with exact functionals in the valley region indicating 
that a precise generalization to two-level models may in 
fact not exist. 

In view of this problem, it is important to realize that 
the main finding of our paper is likely to be insensitive 
to (weak) interaction effects beyond our GGA treatment. 
The reason is, that the statements already follow from a 
two-level model with basic ingredients closed shell system 
(i.e. no magnetism), time-reversal symmetry, inversion 



We have presented a detailed study of charge transport 
properties of the Ceo-molecule coupled to Au-electrodes. 
Our main finding is that the electrical conductance of 
the molecule is strongly suppressed due to two interfering 
transmission channels. The phenomenon was interpreted 
as a precursor to a Fano- anti-resonance. This result has 
been established by combining ab-initio transport cal- 
culations with a toy- model analysis. This analysis also 
suggests, that a thermally driven current is significantly 
less sensitive to such interference effects due to cancella- 
tion effects in the transmission function and the Seebeck 
coefficient. 
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Appendix A: Coefficients for the Grimme correction 



We e mplo y the Grimme empirical correction (Eq. |Al[ 
A2 and A3) 50 to the total GGA energy in order to take 



van-der-Waals interactions into account. The coefficients 



that we used are given in Table III For C and Cu- atoms, 
the values for Cq and Rq have been taken from Ref. [50l 
For gold atoms, the Ro coefficient has been obtained 
from the radius of the electron density contour of a single 
gold atom and the Cq coefficient has been obtained from 
a fit to data obtained from second order Moller-Plesset 
perturbation theory (MP2) 51 . Following Grimme^ we 
used d = 20 and sq — 1.05. 
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-Edisp - ~S& ^ ~^~fdmp(Rij) (Al) 

i=l j=i+l n ij 



fdmp(Rij) ~ 1 + e _ d(iiij/Ro _i) 



(A2) 



(A3) 





C 6 [J.nnAmoP 1 ] 


Ro[A] 


c 


1.75 


1.32 


Cu 


10.8 


1.42 


Au 


21.9 


1.58 



TABLE III: Coefficients used for Grimme empirical correc- 
tion. 



Appendix B: Transport code 

A detailed description of our transport simulations has 
been given in Ref. |29j We present a brief summary of the 
main steps. The effective Kohn-Sham (KS) Hamiltonian 
Hks of the extended molecule is constructed from the 
KS-orbitals and energies previously computed by DFT 
(employing the TURBOMOLE package 34 in our case). 
The Green's function of the extended molecule is built 
from this Hamiltonian i^xs 



G(E) 



E — H K s — X 



where E is the self-energy of the reservoirs, 
taken in the forrrP^ 

s \Snm(Se-iri) m,n e S 

[0 m, n ^ S. 



(Bl) 



It can be 



(B2) 



This self-energy is diagonal in the atomic basis and has 
non-zero values only in the subspace S associated to the 
outermost atomic layer of each electrode. The parame- 
ter r] = 2.72 eV is the leakage parameter. It is adjusted 
such a way that the transmission functions are (approxi- 
mately) invariant under variation of rj. The energy shift 
Se is tuned for each calculation so that the Fermi level 
for the entire system remains unchanged compared to 
what has been obtained from DFT calculation, i.e. about 
the Fermi energy of the metal used for the electrodes. 
Here, we used — 1.25eV < Se < — 1.15eV for Au and 
-2.00eV < Se < -1.94eV for Cu. 

The local density of states (LDOS) projected on the 
molecule can then be written 



LDOS(£) 



(B3) 



where M is the subspace associated to the atoms of the 
molecule. 

The transmission function is given by the following ver- 
sion of the Landauer formula 52 



T(E) — Tr {TlGTrG^} 
with Tl,r defined by 

r L , R = i(s L ,R - 4, R )- 



(B4) 
(B5) 



Appendix C: LDOS projected on molecular orbitals 

In order to illustrate more clearly the mixing of molec- 
ular states, we here introduce a projection of the LDOS 
on specific "coupled molecular orbitals" , The formal 
definition of such orbitals is given by the following con- 
struction. Again, we construct the KS-Hamiltonian Hks 
from the KS-orbitals and energies of the DFT-calculation 
that has been done for the extended molecule (Ceo phis 
parts of the leads). In the orthogonalized atomic basis 
set one can partition Hks m the following way, 



Hks 



Vlc ^rc 

LR ^RC 



(CI) 



where as usual L, R refer to the Hilbert spaces of the left 
and right electrode and C comprises the remaining part 
of the full Hilbert space that belongs to the molecule. 
The states n) are the eigenvectors of the central block 
Hq- They are related to the molecular states of Ceo in 
vacuum, but some effects of the electrode coupling are 
taken into account. Since the states form a complete 
basis of the molecular Hilbert-subspace corresponding to 
Hq, the LDOS on the molecule can then be decomposed 
into the contributions of each molecular orbital 



nr f * 



(C2) 



neM 



The LDOS projected on a molecular orbital can then 
be identified as — G(E)\/j). 

The local density of states projected on the two rele- 
vant molecular orbitals is depicted in Fig. [8] It shows 
that in the valley region of the transmission between 
HOMO- and LUMO-resonances two associated orbitals 
contribute similarly to the LDOS. This suggests that 
there is a possibility for these orbitals to give interfer- 
ing terms in the transmission function T(E). 



Appendix D: Two level (toy) model 

We recall properties of the two-state (toy) model 
that accounts for the transport characteristics of non- 
interacting quantum dots with two effective transport 
levels. As opposed to earlier work 53 , we investigate the 
model analytically in its full parameter space. 
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FIG. 8: (Color) Local density of states projected on the two 
strongly coupled molecular orbitals (MO 1 and MO 2) for the 
geometry where the adatom sits on-top of a C-atom. 




FIG. 9: Schematic representation of the two- level model with 
the inner part representing the Hamiltonian H (grey shaded) 
and the couplings to the leads. 



1. Definition 



Quite generally, the tran smission is given by a formula 
of the Landauer type (|B4|, T(E) = Tr {r L Gr R Gt}. For 
toy model IY,r and G are 2 x 2-matrices. We have for 
the inverse Green's function 



G- 1 (E) = E-H-Z(E) 



(Dl) 



where H is a non-interacting two-level Hamiltonian and 
T,(E) denotes the self-energy that describes the coupling 
to a left and right single chanel wire: E = Er + El- This 



self energy has the general structure 



EaCE) = 9a(E) 



t* 

UOL 
t* 



t* t 




(D2) 



9a(E) (t ua ,tda), a=L,R(D3) 



where t ucn td a denote the hybridization matrix elements 
that connect the two level, up and down, with the wire 
reservoirs; g a (E) resembles a scalar, complex valued 
function, the "surface Greensfunction" of each wire at 
the point contacting to the two-level system. We have 



r a = -i(E a (£7)-Et (£?)); 



(D4) 



Interference effects can occur, if I\r and G do not com- 
mute, so that they cannot be diagonalized simultane- 
ously. In order to highlight them, we begin substituting 
( D3 ) into the trace formula; we obtain 



T(E) = (27r) 2 g L g R 



(tuLtdh)G 



l uR 

t* 
L dR 



(D5) 



with the contact density of states g(E) = -^$sg a (E). 
Next, we rotate into the basis of eigenfunctions of G{E) 



G(E) = U I E ~ z ° 





E-z lt 



(D6) 



In general, the pole positions zq,i an d the eigenvectors 
comprising the columns of the rotation matrix U inherit a 
dependency on energy, E, through E(E'). In order to sim- 
plify notation we introduce effective hybridization matrix 
elements, v 



U 0R 



J 1R 



and transmission coefficients 

Ta(3(E) = 27ry/g a g^(t ua tda)G 

E — zq E — z\ 
that allow us to write 

T(E) = \r LR (E)\ 2 



t* 

l uR 

t* 
L dR 



(D7) 
(D8) 



L u(3 
l d(3 



(D9) 



(D10) 



To explicitly single out the interference term, we employ 
a decomposition 



T(E) = T (E) + T^E) + T m {E). 



(Dll) 



The two first terms constitute the non-mixing contribu- 
tions from each energy level 

rp ( ^ |^ol|>or| 2 Kl| 2 Kr| 2 m , 

To{E)= \E-z \i ' Tl{E)= \E- Zl \* (D12) 
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Each term refers to a single pole only and thus is re- 
produced by the model of isolated resonances, Eq. |2j. 
Interference enters via the mixed term 



Toi(E) 
G(E) 



2^/^008 (* c -e) 

arg((£ - z )(E - z{)) 



(D13) 
(D14) 
(D15) 



As is seen from equation (D13), interference effects are 
controlled by two angles, \I/ C and O. They exhibit quite 
different generic properties. \£ c carries an energy depen- 
dency only via T>(E) that reflects a dispersion in the (lo- 
cal) density of states in the reservoirs. Because the lat- 
ter often is very smooth compared to the level splitting, 
|eo — ei| with $lzi = €i, it can typically be ignored, so that 
for practical purposes \£ c is energy independent. By con- 
trast, Q(E) can exhibit a very sharp energy dependency, 
especially in the limit of low damping. 



Remark: Unitarity theorem 

We split the transmission coeffients into hermitian and 
anti- hermit ian contributions: 



Tap = 27T yj Q a Qp (t ua tda ) G ™ 



u(3 



l dj3 



i(G+Gt)+^GrGt 



Z up 

t* 
b dp 



The first term simplifies to \ 
ond one takes the form 



Ta[3 +T (3a 



while the sec- 



(27T) 2 



2i 



\fQaQp{t ua tda) 



■ G 



- 1 V" 



{tuatda) 



da 



q] I tup 



(D23) 



Remark: Level broadening 

To reveal further transport properties of the two-level 
transmission function, we relate the amplitudes \vo a \ 2 to 
the level broadenings 7^. This broadening is originates 
from the anti-hermitian part of the inverse Green's func- 
tion 



G" 1 = Go 1 + |r (D16) 
Go 1 = E-F-i(E + St). (D17) 



We have 



TrG -1 = 2E - z - z x = Tr Gq 1 + -Tr T. (D18) 



Since G 1 , T by construction are hermitian, each trace is 
real; hence 



K[2b + «i] = Tr 
Z[z + Zl } = ^TrT. 



fr + -(E + st) 



(D19) 
(D20) 



The second line simplifies after recalling (D3): 



%[z + z 1 ] = ^-2^^p a (|^| 2 + |^| 2 )^(D21) 

= -J^K*| 2 + K*I 2 



(D22) 



Summarizing, we have in a matrix notation the general 
statement 



T _ r t — _ \ TT \ 



(D24) 



that satisfies a unitarity theorem r 1 — [r^] 1 = i. 



2. Time reversal symmetry 

In the presence of time reversal symmetry the matrices 
Gq 1 and T are both real symmetric, so G _1 is (complex) 
symmetric and U is (complex) orthogonal. Then simpli- 
fications arise for Tbi. 

We focus the discussion on the weak coupling limit, 
where the resonance positions are split by an amount that 
considerably exceeds their broadening, |ei— eo| ^> 7o,i- In 
this case we can consider the anti-hermitian piece of the 
self energy as a perturbation and the eigenvectors 

are real to first order in 7^ at least to this order U is 
real orthogonal. Hence, also the effective hybridization 
matrix elements Vi a are real and therefore \I/ C takes values 
zero or 7r; we have 

cos(^ c -G(E)) = sign(v 0L v 0R v 1L v m ) cosG(E) (D25) 

To the leading order in 7^ the energy dependency is de- 
scribed by 



cos 0(E) = sign ([E-eoHE-ei]). 



(D26) 



Equations ( D25|D25| ) are important because they high- 
light two characteristic features of the interference term 
Tqi. First, whether constructive (cos(\I/ c — 6) > 0) or 
destructive (cos(\I/ c — 6) < 0) interference prevails in 
the intermediate range €q < E < e± is controlled by the 
model dependent first factor in ( |D25 ). If it is positive, 
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\I/ C = 0, interference is destructive; it is constructive in 
the alternative case \I/ C = tt. 

Second, Tqi has a very sharp dependency on energy 
when E sweeps by eo or e\. In fact, the corresponding 
deriv ative diverges in the limit ji —> 0; specifically, with 
Eq. (D15) we have 



cos 0(E) 
sin6(£;) 



(E - e )(E - ei) +7o7i 
\E-Zol\E-zt\ 

7i(g ~ eg) + 7o(^ ~ gi) 
LE-*b|LE-*i| 



(D27) 
(D28) 



3. 



Limiting cases: symmetric vs asymmetric lead 
couplings 



We now analyze two limiting cases in which the angle 
\I/ C is either close to tt (constructive interference) or close 
to zero (destructive interference). 



energy dependency of the eigenvalues is weak (i.e. in the 
wide band limit, where g a (E) is nearly constant), Eq. 
(D30) has the real numbered root 



E* 



+ h\v \ 2 



\vo\ 2 



N 2 



(D31) 



It constitutes a weighed average that is situated between 
the two pole positions. 

Since at E* we have complete destructive interference, 
it is clear that at this energy \I/ C — Q(E*) = tt. Since in 
the valley region is close to tt we have <£ c « with 
corrections in 70,1. Therefore, in the valley region away 
from the anti-resonance E* we expect 



Tqi =2 v ^T 1 cosO(E) 



(D32) 



with cosO(i^) as given in Eq. (D27). 



a. Fully symmetric coupling: Fano anti-resonance 

We begin with the second situation, destructive inter- 
ference, which is the easier one to investigate and exploit 
Eq. (D24) now for a system with completely symmetric 
coupling, t uL = t uR = t u , t dL = t dR = t d . This situation 
is realized, e.g., to a very good approximation for cross- 
conjugated molecular wires with side coupling chains as 
investigated in Ref. [54 . A cross-conjugated molecule 
was studied recently experimentally in Ref. [55], where 
the Fano- ant iresonance has indeed been observed. 

Our interest is in the off-diagonal matrix element, 

tlr(£) = K[r LR (^)]-^^±^|r LR (E)| 2 {D29) 

This equation implies the following fact: Let E* be a 
root of the real part of the transmission coefficients, 
3?[^lr(^*)] = 0. Then, at this energy also the imaginary 
piece has a vanishing physical solution: ^[tlr(£^*)] = 0. 

For generic situations E* can be shown to be real, so 
that there is no transmission at this energy, T(E*) = 0. 
Namely, let U(E) be the unitary rotation diagonalizing 
the hermitian matrix G + ; then 



5K[tlr(£)] = 2^ v /^ 



M 2 



N 2 



E - e E-ei 



(D30) 



where eo,i denote the real eigenvalues of (G + G^)/2 and 
(vq,vi) = (t u ,td)U(E). Under the assumption, that the 



b. Fully asymetric coupling 

The previous example suggests, that constructive in- 
terference could be expected for an asymetri c limit , 
t U R,tdL = 0. In this case, because of equations ( |D7|D8 ) 
the effective hybridization matrix elements reproduce the 
entries of U and Z7 -1 : 



J iR 



so that 



* c = arg (C/oo^-'loi^-'ln^oi, 



(D33) 



The product in brackets is easy to evaluate in the pres- 
ence of time reversal symmetry recalling that U is or- 
thogonal in this case. Like any orthogonal 2 x 2-rotation 
matrix it has a representation 



U 



cos w sin w 
— sin w cos w 



w e C 



We obtain \I/ C = n. Hence, we now expect a constructive 
interference contribution 



Tqi = -2^y^V^ 1 cosG(E). 



(D34) 



where again cos0(^) is given in Eq. (D27). 



X. Lu, M. Grobis, K. H. Khoo, S. G. Louie, and M. F. 



Crommie, Phys. Rev. Lett. 90, 096802 (2003), URL http: 



//link . aps . org/doi/10 . 1 103/PhysRevLett . 90 . 096802 



2 X. Lu, M. Grobis, K. H. Khoo, S. G. Louie, and M. F. 
Crommie, Phys. Rev. B 70, 115418 (2004). 

3 C. Rogero, J. I. Pascual, J. Gomez-Herrero, and A. M. 



13 



Baro, J. Chem. Phys. 116, 832 (2002). 
4 J. D. Sau, J. B. Neaton, H. J. Choi, S. G. Louie, and M. L. 



Cohen, Phys. Rev. Lett. 101, 026804 (2008) , URL|http:| 
//link . aps . org/doi/10 . 1103/PhysRevLett . 101 . 02 6804 
° I. F. Torrente, K. J. Franke, and J. I. Pascual, J. Phys. : 
Condens. Matter 20, 184001 (2008). 

6 A. Tamai, A. P. Seitsonen, F. Baumberger, M. Hengs- 
berger, Z.-X. Shen, T. Greber, and J. Osterwalder, Phys. 
Rev. B 77, 075134 (2008). 

7 E. Abad, C. Ganzalez, J. Ortega, and F. Flores, Organic 
Electronics 11, 332 (2010). 

8 I. Hamada and M . Tsukada, Phys. Rev. B 83 
245437 (2011), URL |http: //link. aps.org/doi/ 10. 1103/| 
PhysRevB . 83 . 245437| 

y H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. 
Alivisatos, and P. L. McEuen, Nature 407, 57 (2000). 

10 L. H. Yu and D. Natelson, Nano Letters 4, 79 (2004). 

11 A. N. Pasupathy, R. C. Bialczak, J. Martinek, J. E. Grose, 
L. A. K. Donev, P. L. McEuen, and D. C. Ralph, Science 
306, 86 (2004). 

12 J. J. Parks, A. R. Champagne, G. R. Hutchison, S. Flores- 
Torres, H. D. Abruna, and D. C. Ralph, Phys. Rev. Lett. 
99, 026601 (2007). 

13 N. Roch, S. Florens, V. Bouchiat, W. Wernsdorfer, and 
F. Balestro, Nature 453, 633 (2008). 

14 G. Schulze, K. J. Franke, and J. I. P. and, New J. of Phys. 
10, 065005 (2008). 

15 G. Schulze, K. J. Franke, A. Gagliardi, G. Romano, C. S. 
Lin, A. L. Rosa, T. A. Niehaus, and T. Frauenheim, Phys. 
Rev. Lett. 100, 136801 (2008). 

16 T. Frederiksen, K. J. Franke, A. Arnau, G. Schulze, J. I. 
Pascual, and N. Lorente, Phys. Rev. B 78, 233401 (2008). 

17 A. Gagliardi, G. Romano, A. Pecchia, A. D. Carlo, 
T. Frauenheim, and T. A. Nienhaus, New Journal of 
Physics 10, 065020 (2008). 

18 J. Fock, J. K. Sorensen, E. Lortscher, T. Vosch, C. A. 
Martin, H. Riel, K. Kilsa, T. Bjornholm, and H. van der 
Zant, Phys. Chem. Chem. Phys. 13, 143 25 (2011), URL 
|http : //dx . doi . org/10 . 1039/C1CP20861F1 

19 N. Neel, J. Kroger, L. Limot, T. Frederiksen, M. Brand- 
byge, and R. Berndt, Phys. Rev. Lett. 98, 065502 (2007). 

20 G. Schull, T. Frederiksen, M. Brandbyge, and R. Berndt, 
Phys. Rev. Lett. 103, 206803 (2009). 

21 G. Schull, T. Frederiksen, A. Arnau, D. Sanchez-Portal, 
and R. Berndt, Nature Nanotechnology 6, 23 (2010). 

22 G. Schull, Y. J. Dappe, C. Gonzalez, H. Bulou, and 
R. Berndt, Nano Letters 11, 3142 (2011). 

23 T. Bohler, A. Edtba uer, and E. Scheer, Phy s. Rev. B 76, 
125432 (2007), U RL|http : //link . aps . org/doi/10 . 1103/| 
PhysRevB. 76. 125432| 

M. Kiguchi and K. Murakoshi, J. Phys. Chem. C 112, 8140 
(2008). 

25 S. K. Yee, J. A. Malen, A. Majumdar, and R. A. Segalman, 
Nano Letters 11, 4089 (2011). 

26 M. Mayor, H. B. Weber, J. Reichert, M. Elbing, C. von 
Hanisch, D. Beckmann, and M. Fischer, Angwandte 
Chemie-International Edition 2003, 5834 (2003). 

27 C. A. Martin, D. Ding, J. K. Sorensen, T. Bjornholm, 
J. van Ruitenbeek, and H. S. J. van der Zant, J. Am. Chem. 
Soc. 130, 13198 (2008). 

28 E. Leary, M. T. Gonzalez, C. van der Pol, M. R. 
Bryce, S. Filippone, N. Martin, G. Rubio-Bollinger, and 
N. Agrait, Nano Letters 11, 2236 (2011). 

29 A. Arnold, F. Weigend, and F. Evers, Journal of Chemical 



Physics 126, 174101 (2007). 

M. K. Shukla, M. Dubey, and Leszczynski, ACS Nano 2, 
227 (2008). 

M. K. Shukla, M. Dubey, E. Zakar, and J. Leszczynksi, J. 
Phys. Chem. 113, 11351 (2009). 

S. Bilan, L. A. Zotti, F. Pauly, and J. C. Cuevas, 
arXiv:1203.3101vl (2012). 

T. Ono and K. Hirose, Phys. Rev. Lett. 98, 026804 (2007). 

R. Ahlrichs, M. Bar, M. Haser, H. Horn, and C. Kolmel, 

Chem. Phys. Lett. 162, 165 (1989). 

A. Becke, Phys. Rev. A 38, 3098 (1988). 

S. Grimme, Journal of Computational Chemistry 25, 1463 

(2004). 

S. T. Schneebeli, M. Kamenetska, Z. Cheng, R. Skouta, 
R. A. Friesner, L. Venkataraman, and R. Breslow, J. Am. 
Soc. Chem. 133, 2136 (2011). 

R. W. Lof, M. A. van Veenendaal, B. Koopmans, H. T. 
Jonkman, and G. A. Sawatzky, Phys. Rev. Lett. 68, 3924 
(1992). 

U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986), URL 
|http : //link . aps . org/doi/10 . 1 103 /PhysRevB . 33 . 551 
S. H. Ke, M. Yang, S. Cortarolo, and H. U. Baranger, Nano 
Lett. 9, 1011 (2009). 

J. R. Widawsky, P. Darancet, J. B. Neaton, and 
L. Venkataraman, Nano Lett. 354, 12 (2012). 
M. Pau lsson and S. Datta, Phys. Rev. B 67, 241403 (2003), 
URL |http : //link . aps . org/doi/10 . 1 103/PhysRevB . 67 . | 
24140 



D. Segal, Phys. Rev. B 72, 165426 (2005), URL|http:| 
//link. aps.org/doi/10.1103/PhysRevB. 72. 165426| " 
^ S.-H. Ke, W. Yang, and H. U. Baranger, Nano Lett. 8, 
3257 (2008). 

45 D. M. Cardamone, C. A. Stafford, and S. Mazumbdar, 
Nano Lett. 6, 2423 (2006). 

46 J. P. Bergfield, Z.-F. Liu, K. Burke, and C. A. Stafford, 
Phys. Rev. Lett. 108, 066801 (2012). 

47 G. Stefanucci and S. Kurth, Phys. Rev. Lett. 107, 216401 
(2011). 

48 P. Schmitteckert and F. Evers, Phys. Chem. Chem. Phys. 
13, 14417 (2011). 

49 P. Schmitteckert and F. Evers, Phys. Rev. Lett. 100, 
086401 (2008). 

50 S. Grimme, Journal of Computational Chemistry 27, 1787 
(2006). 

51 M. Peter, Master's thesis, Universitat Karlsruhe (2009). 

52 Y. Meir and N. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

53 G. C. Solomon, D. Q. Andrews, T. Hansen, R. H. Gold- 
smith, M. R. Wasielewski, R. P. Van Duyne, and M. A. 
Ratner, J. Chem. Phys. 129, 054701 (2008). 

54 G. C. Solomon, D. Q. Andrews, R. H. Goldsmith, 
T. Hansen, M. R. Wasielewski, R. P. Van Duyne, and M. A. 
Ratner, J. Chem. Phys. 130, 17301 (2008). 

55 C. M. Guedon, H. Valkenier, T. Markussen, K. S. Thyge- 
sen, J. C. Hummelen, and S. J. van der Molen, Nat. Nano. 
7, 305 (2012). 

56 Here, we use a slightly simplified version in which we 
choose an approximation ^ c = 0, 7r. In this way we elimi- 
nate two fitting parameters, the real angle ^ c and a back- 
ground contribution that cuts off the Fano- anti-resonance. 
Both parameters need to be accounted for if the simplified 
(toy) model, when adjusted near the resonance positions, 
does not capture the behavior of the transmission function 
in the valley region. 



